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^ ■ Abstract 

. The absorption of a high-energy photon from the external cosmic gamma-ray background in 

Q \ the inner neutron star magnetosphere triggers the generation of a secondary electron-positron 

H ■ 

' plasma and gives rise to a lightning — a lengthening and simultaneously expanding plasma tube. It 

: 

propagates along magnetic fields lines with a velocity close to the speed of light. The high electron- 

\ positron plasma generation rate leads to dynamical screening of the longitudinal electric field that 

Q>^ I is provided not by charge separation but by electric current growth in the lightning. The lightning 

■ radius is comparable to the polar cap radius of a radio pulsar. The number of electron-positron 

. pairs produced in the lightning in its lifetime reaches 10^*^. The density of the forming plasma is 

\ comparable to or even higher than that in the polar cap regions of ordinary pulsars. This suggests 

^ ' that the radio emission from individual lightnings can be observed. Since the formation time of 

ixl ■ 

sj . the radio emission is limited by the lightning lifetime, the possible single short radio bursts may 
^\ 

be associated with rotating radio transients (RRATs). 
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1. INTRODUCTION 



This paper is a continuation of our previous paper [1], in which we calculated the source 
of electron-positron pairs for the case of a strong electric field in a vacuum neutron star 
magnetosphere. Interest in this subject matter is related to the appearance of new obser- 
vational data on nonstationary radio sources, such as intermittent pulsars [2] and rotating 
radio transients (RRATs) [3, 4] . These sources, just as ordinary radio pulsars, are associated 
with neutron stars. In this case, the corresponding surface magnetic field is high and reaches 
5 X 10^^ G [3, 5], i.e., it exceeds the Schwinger field, for example, for RRAT J1819-1458. 
Magnetars, which are also characterized by radio emission nonstationarity observed at low 
[6-9] and high [10-12] radio frequencies, have an even stronger magnetic field. Note that 
the presence of such a strong magnetic field and an electron-positron plasma changes the 
kinematics of quantum processes [13, 14]. For example, photon splitting becomes significant 
[15, 16]. However, it does not suppress the generation of an electron-positron plasma in 
a superstrong magnetic field and, consequently, the corresponding radio emission [17, 18], 
although the radio emission itself can be observed erratically. In particular, there are data 
on the switch-off of a magnetar without any accompanying X-ray variability [12]. Since 
plasma outflows from the magnetosphere are currently believed to be responsible for the 
observed radio emission from neutron stars (see, e.g., the review [19]), the activity of all 
the mentioned radio sources suggests a possible cessation of the plasma generation in their 
magnetospheres [20]. This is also backed by the measured difference in spin-down rate of 
the intermittent pulsars PSR B1931+24 and PSR J1832+0029 in the periods of "operation" 
and "silence" [21]. If, alternatively, there is no plasma and no free charge escape from the 
stellar surface occurs, which is confirmed by the spin-down in the period of pulsar silence 
[20], then the magnetosphere can be filled through the absorption of photons from the ex- 
ternal cosmic gamma-ray background [22-24] . Allowance for the influence of these photons 
on the particle production processes was shown to be necessary in [25]. It was pointed out 
in the latter paper that an influx of new charges that can be provided by the absorption 
of photons from the diffuse gamma-ray background is needed in the Ruderman-Sutherland 
model [26] to trigger each new spark in the polar gap of a pulsar magnetosphere. 

Let there be a primary galactic photon whose energy and propagation direction are such 
that the transverse momentum component exceeds 2meC, where mg is the electron mass and 
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c is the speed of light. This photon can then produce an electron-positron pair in a magnetic 
field [27]. The pair particles emit synchrotron photons as a result of the transition to the 
zeroth Landau level and, accelerating in a strong electric field present in the magnetosphere, 
begin to emit curvature photons. The absorption of these photons in a magnetic field leads 
to the production of the next generation of particles. The developing cascade gives rise to a 
lightning — a plasma tube propagating almost along magnetic field lines. In this paper, we 
study in detail the formation of this tube in a neutron star magnetosphere and the generation 
of an electron-positron plasma in it. 

The paper is structured as follows. In Section 2, we provide the basic equations describing 
the motion of particles and the evolution of the electromagnetic field and pass to complex 
variables in the Maxwell equations. In Section 3, we investigate the dynamical screening of 
the electric field in a lightning. In Section 4, we calculate the self-consistent Lorentz factor 
of particles and the longitudinal electric field at the lightning center. In Sections 5 and 6, 
we calculate the distributions of electrons and positrons and the total number of produced 
electron-positron pairs, respectively, before the onset of screening and during the electric 
field screening. In Section 7, the calculation is concretized for a power-law dependence of 
the electric field on the radial coordinate. In Section 8, we estimate the radius of the forming 
lightning and determine the multiplicity parameter. Our main results and conclusions are 
summarized in Conclusions. 



2. THE MOTION OF THE PRODUCED PARTICLES AND THEIR INFLUENCE 
ON THE ELECTROMAGNETIC FIELD 

The motion of electrons and positrons in a neutron star magnetosphere is ultrarelativistic. 
Because of the presence of a strong accelerating electric field, the electrons and positrons 
move in opposite directions. We can write the relation 

v± = ±b-Fexb, (1) 

where v+ is the positron velocity, v_ is the electron velocity, b = B/i? is a unit vector 
directed along a magnetic field line, and e — Fi/B is the electric field vector normahzed to 
the magnetic field strength. Formula (1) describes the ultrarelativistic motion of charged 
particles along magnetic field lines and the drift motion across them. The described consid- 
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eration is proper even in Deutsch's unscreened vacuum electromagnetic field [28] , while the 
centrifugal drift may be neglected [22]. 

The continuity equations for the electron and positron densities are 

^ + V • n± v± = g(n+ + n_) , (2) 

where n+ and n_ are the positron and electron number densities, respectively, Q is the total 
number of electron-positron pairs produced per unit time per particle. These are the basic 
equations for describing the generation of an electron-positron plasma. 

Here, as in our previous paper [1], wc will use a dimcnsionless system of units. We will 
measure the electric and magnetic field strengths in units of the critical field 

B^^ = !^ ^ 4AU X W G, 
eh 

where e is the positron charge. The length and time will be measured in units of the 
Compton electron wavelength X = h/nieC « 3.862 x lO"-*^^ cm and its ratio to the speed of 
light X/c^ 1.288 X 10~^^ s, respectively. Note the useful relations 1 cm pa 2.590 x 10^° and 
1 s 7.763 X 10^0. 

The Maxwell equations in dimensionless units in cylindrical coordinates r, (p, z are 



1 d(rEr) BE, ^ 

+ o = 47rQ;pe, 
r or oz 


(3) 


dE^ dBr 
dz dt ' 


(4) 


dEr dE, dB^ 
dz dr dt ^ 


(5) 


1 d{rE^) _ dB, 

r dr dt ^ 


(6) 


ld{rBr) 1 dB, _^ 
r dr dz ' 


(7) 


dB^ ^ . , dEr 
dz ='""^^+ dt ' 


(8) 


dBr dB, 1 dE^ 
dz dr dt ^ 


(9) 


r dr dt 


(10) 



Here, z is the longitudinal coordinate measured along the magnetic field direction, r is 
the radial coordinate, </? is the azimuth angle, t is the time, E and B are the electric and 
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magnetic field strengths, pe — n+ — n_ and j = n-+v+ — n_v_ are the charge and current 
densities, and a — /he ~ 1/137 is the fine-structure constant. In view of the assumed 
local cylindrical symmetry of the plasma tube, we set d/dip = 0. Since the forming plasma 
tube will be curved due to the existing curvature of the magnetic field lines, these equations 
are applicable under the condition R <ti p, where R is the plasma tube radius and p is the 
radius of curvature of the magnetic field lines. At the end of the paper, we show that the 
tube radius is limited by 100 m (see Eq. (68)), while the radius of curvature of the magnetic 
field lines exceeds the neutron star radius Rs ~ 10 km. Consequently, this condition is 
always met. 

Let us calculate the electron and positron velocities. For this purpose, let us turn to 
Eq. (1). The particle velocity is the sum of the longitudinal component giving the motion 
along a magnetic field line and the transverse component determined by the electric drift in 
crossed fields. In order of magnitude, the former and latter components are equal to 1 and 
E/B, respectively. Initially, as long as the number of produced particles is small, the total 
electromagnetic field is specified by the external electric field E^^* and the external magnetic 
field B^^*. When, however, the plasma tube formation begins, the internal electromagnetic 
field also emerges. Let us specify the internal electric field E*"* and the internal magnetic 
field B"**. All scalar components of these fields are equal between themselves in order of 
magnitude. The forming internal electric field will screen the external longitudinal electric 
field, but, obviously, it will be unable to exceed it in absolute value. We can write S'"* ~ 
^mt < ^ext Qj.^ passing to the components of the total field that determines the particle 
velocities, Br B^ ^ Er ^ ^ ^ Ef^* and Bz ~ Bf^*. In addition, we should 
take into account the fact that the longitudinal component Bz of the total magnetic field 
prevails over the orthogonal components B^. and B^p. For this reason, the total magnetic 
field B coincides with the component B^ to within terms quadratic in B^/Bz ~ B^/Bz- 
Using these estimates, expanding the expressions in Eq. (1) in terms of E^^^/Bf^* to the 
first order inclusive, and discarding the terms of the next orders, we will write the electron 
and positron velocity components 

= ±br + e^, = ±b^ - Cr, vf = ±1. (11) 

The particle velocities and number densities and the source of particles which depends 
on the longitudinal components of the electric and magnetic fields, enter into Eq. (2). Apart 
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from the longitudinal components, the orthogonal components of the fields, which are not 
required per se, enter into the Maxwell equations. To exclude them from consideration, let 
us introduce four new sought-for functions 



proportional to the orthogonal particle velocity components: — Bvf and — Bv^. It 
is also convenient to pass from the variables z, t to new variables 

We can now rewrite the Maxwell equations (3)-(10) by specifying the functions V^, V~ , 
V^, V~ , Ezi and B^ dependent on the variables i+,i_,r as unknown quantities: 

dV^ . dB 



2^ = -4™j,-^, (12) 

1^ . (13) 

r or ot± 

dV^ dE 

l^ = ±4.«(,.„J±2fi. (15) 

r or at± 

Let us introduce a complex electromagnetic field vector F = E+iB, a complex transverse 
velocity = + iv^, and a complex transverse electric current density j± = jr + ij,p- 
Equations (12)-(15) will then take a more compact form: 



2^ = .^4™j,± ^ 



1 d{rVl 



(QP± \ 
47ran^ + _^j, (16) 



r dr 

where we introduced the quantity = Bv"^ and the longitudinal component = F^ of the 
complex electromagnetic field vector and the corresponding complex conjugate component 
F~ — F*. It is also easy to derive the equations for the longitudinal components F^ that 
we will need below: 



(17) 



They are convenient in that they contain the dependence of the longitudinal component of 
the complex electromagnetic field only on the electron and positron densities and the complex 
transverse current density. However, the transverse current is not specified independently 
but is expressed in terms of the particle densities and velocities as follows: 

h^^{n+Vl-n^Vz). (18) 

Below, we will expand all of the quantities entering into the equations in terms of the 
radial coordinate r. If we arbitrarily designate some sought-for quantity as y, then its 
expansion in general form is represented by a power series 



where the quantities 



n=0 



y 



(n)_ ^ 



r=0 

depend only on the coordinates i or t- and do not depend on r. Obviously, the positron 
and electron densities, the charge and current densities, and the longitudinal electromagnetic 
field components are expanded in terms of even powers of r, while the transverse components 
of the electromagnetic field and the particle velocities are expanded in terms of odd powers 
of r. For example, 

= FW + Fi2)^ + 0(r^), n± = nf +0(r2), v^ = vf'-\ ^ 0{r'). 

Expanding all of the quantities entering into the equation into series of r and equating the 
terms with the same powers of r, we will obtain an infinite set of equations among which 
we will restrict our consideration to the lowest-order equation. 

To find the expansion of the current density we will need the expansion of Eq. (16): 

±(1) _ I A^^,A0) 9F^^^^ 



Now, using Eq. (18), we easily obtain 

.(1) i ( (o)dFl^ {o)dFV^ (o) (o) 

= l^n^ ^ + n_ + 8™+ n_ 

The continuity equations for finding the positron and electron densities 

2^ + i^^^ = Qj. (19) 

at± r or 



after the expansion take the form 

^ + nL^.±W = ig(°)jf. (20) 

The equation for the longitudinal electromagnetic field component follows from the expan- 
sion of Eq. (17): 

2Fi^) + □,, = S.a(^^ + ^- j , (21) 
where we introduced the D'Alembert operator in longitudinal coordinate 



Now, it remains to express the derivatives of the positron and electron densities in terms 
of the particle source and the radial fluxes using Eqs. (20) and to substitute the derived 
expressions into Eq. (21). Separating out the imaginary part of the complex current density 
j_L allows the azimuthal current density to be obtained: 

.(1) 1 / (G)dEf'' (o)dEi°^ (0) (0) 

The expansions of the radial positron and electron velocities can be derived from the relation 
v^^^^ — V^^^^ / B^'^\ once the real part of V^^^^ has been separated out: 

y±(i) ^ 

Finally, we obtain the equation for the longitudinal component of the complex electromag- 
netic field vector 



2F(2) + nilFi°) = 87ro; 



^(0) -(0) / (o)<9Fi°^ {o)dF^°^ (0) (0) 



(22) 



This equation is convenient in that it contains only the longitudinal components E^^^ and 
of the electromagnetic field and the electron and positron densities. The radial and 
azimuthal components of the electromagnetic field and the particle velocities turned out to 
be excluded. 



3. DYNAMICAL SCREENING OF THE LONGITUDINAL ELECTRIC FIELD 

It was shown in [29] that the first term on the right-hand side of Eq. (22) prevails over 
the remaining terms. Note that this can be verified only after the calculation of the self- 
consistent source of electron-positron pairs by taking into account the dynamical screening 
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of the external longitudinal electric field. Let us turn to the left-hand side of Eq. (22). 
Initially, this expression is equal to OFz. The D'Alembert operator is represented as the 
sum 

□ = A^ + Dii, 

where the two-dimensional Laplace operator in radial coordinate was introduced: 

1 d d 

A_L = r — . 

r dr dr 

Let the characteristic longitudinal distance and the characteristic time scale on which the 
electromagnetic field changes be Z and T, respectively. The field changes with radial coor- 
dinate on distances of the order of the radius it! of the forming plasma tube. We can then 
symbolically write 

1 ^11 

/\ I II — 

^ i?2' II ^2 2-2- 

We will assume that, first, the change of the electromagnetic field is quasi-stationary and, 
second, the characteristic longitudinal scales of the field change exceed considerably the 
corresponding transverse scales. Then, this means that Dy <^ Aj^. The stipulated condi- 
tions allow us to formally set Dn = in Eq. (22) and to obtain the final equation for the 
longitudinal component of the complex electromagnetic field: 



Once the real and imaginary parts have been separated out, we have 

=47raQ(°)jf, = 0. (23) 

The radial distribution of the longitudinal electric field inside the plasma tube is 



^2 



Let us specify some electric field E^"^"-'^ = Ez\r=R at the tube boundary. To determine this 
field, we should take into account the fact that the longitudinal electric field at infinity is 
equal to the external one, Ez\r=oo — Next, let us turn to Eq. (5). The radial derivative 
dEz/dr entering into it exceeds the derivatives dEr/dz and dB^/dt in order of magnitude. 
This follows from the field estimates pointed out after Eq. (10) and from the inequahties 
Z ^ R and R <^T. We arrive at the equality dE^/dr — whose integration leads us to 



the boundary condition E^'^'"^ — Ef^^. Thus, the longitudinal electric field at the plasma 
tube boundary is approximately equal to the external longitudinal electric field existing in 
the magnetosphere in the absence of an electron-positron plasma. The longitudinal electric 
field at the tube center is then 

Ef) ^ET' -Ef^ — . (24) 
2 

Similar reasoning is also valid for the longitudinal magnetic field. We should only turn 
to Eq. (9) and take into account the fact that j^p = outside the plasma tube. The 
boundary magnetic field B^°^'^^ — Bz\r=R then coincides with the external longitudinal 
magnetic field The equahty = ^bound alternatively, 

^(0) ^ ^ext (25) 

follows from Eq. (23). 

Consider a quasi-stationary plasma tube. Obviously, the longitudinal electric current 
existing in the plasma tube gives rise to a nonzero azimuthal magnetic field B^. Since the 
particles move mainly along magnetic field lines, the appearance of an azimuthal magnetic 
field component leads to the appearance of an azimuthal particle velocity component v^. 
Here, in general, the existence of an additional drift velocity component due to the presence 
of a radial electric field Er should also be taken into account (see Eq. (11)). This is unimpor- 
tant for qualitative considerations, because we can consider the regions in the plasma tube 
where the electron and positron densities are comparable between themselves but fairly high 
per sc. The radial electric field will then be low, while a certain electric current causing the 
magnetic field lines to twist will flow in the tube. An azimuthal electric current component 

will appear, which, in turn, will cause the longitudinal magnetic field to change. This 
change will be nonuniform in radial coordinate. Consequently, the second derivative of the 
magnetic field is still nonzero. However, it is small compared to the second derivative of the 
longitudinal electric field. This is because the production of a secondary electron-positron 
plasma is fairly cflicicnt due to the large Q. Therefore, Eq. (25) should be understood in 
the sense that the change in longitudinal magnetic field through the appearance of an inter- 
nal screening electromagnetic field in the plasma tube is small compared to the change in 
longitudinal electric field caused by the same internal field. 

Now, returning to Eq. (23), we can understand the longitudinal electric field screening 
mechanism. Here, it is important that the electron-positron pair production rate primarily 
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determines the second spatial derivatives of the longitudinal electric field. In other words, 
the curvature of the two-dimensional cylindrically symmetric surface Ez{r) increases with 
plasma generation rate. Since the boundary condition E^'^'^ — El^* is met, the surface sag 
causes the longitudinal electric field at the plasma tube center to decrease. This mechanism 
may be called dynamical screening, because the surface sag depends on the plasma generation 
rate. It differs from the mechanism of ordinary static screening related to charge separation. 
The existence of screening different from static one was shown in [30]. The electric field 
screening in the sparks forming in the polar gap considered in the latter paper was also 
called dynamical one. 

Let us integrate Eq. (19) over the radial coordinate with the weight 7rr from to oo: 

R 



= TT y Qj.rdr. (26) 



dt± 



Here, and N_ denote the linear positron and electron densities, respectively. As above, 
let us expand all quantities in terms of r. Let us also introduce the linear charge density 
Np — — and the electric current Jz = + . We can then write the relations 

= rzf 7^i?^ Np = p<S>\R', J, = Tri?^ (27) 

Equation (26) transforms to 

Q^'^J. .... 

9^ - 

We can also pass from the variables N± to the variables Np and Jz and obtain 

Thus, to find the distribution of electrons and positrons along the longitudinal coordinate 
and to calculate the number of particles produced in the plasma tube lifetime, the source of 
electron-positron pairs per particle (5^°^ should be found. 

4. THE PARTICLE LORENTZ FACTOR 

Formulas (23), (24), and (27) allow the dynamical screening of the longitudinal electric 
field to be taken into account: 

^(0) =£;|^t_2ag(o)j^. 
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As Q^^^ , we should take the effective local source 

Q'" = (31) 

calculated for the plasma tube center [1]. The source Q^-^-^ depends on the logarithmic factor 
^0 ~ 0.1 and the stationary Lorentz factor 

/ o \ 1/4 

70 ^(^^„/) . (32) 

which, in turn, is a function of the longitudinal electric field E\\ = E^. For the tube center, we 
have Ez = Ei^\ The described dependence allows the Lorentz factor 70 and the longitudinal 
electric field £"1°^ at the plasma tube center to be calculated self-consistently. The Lorentz 
factor 7o can be found from the equation 

< + f ^.7o' - iLt = 0, (33) 

where ^ext is defined by Eq. (32), in which we should set = 

Equation (33) is a quartic algebraic equation. It can be brought to the form (for more 
detail, see [31]) 



(34) 



where we can take 



as K and L. Here, Y is one of the solutions of the cubic equation 

+ DiY + L>2 = 0, (35) 

where 

Since Di > 0, this equation has one real and two complex conjugate roots. Let Y mean the 
real root of Eq. (35). Cardano's formula gives 



3 / 
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where the real values should be taken for the left and right cubic roots. Obviously, Y < 0. 

Let us return to Eq. (34). The quadratic equation obtained when the left bracket is equal 
to zero will not be of interest to us, because, as is easy to verify, it has a negative discriminant. 
This corresponds to the existence of two complex conjugate roots. The quadratic equation 
obtained when the right bracket is set equal to zero has two real roots of opposite signs. 
Choosing the positive root, we obtain the required expression for the Lorentz factor: 



The corresponding total longitudinal electric field at the plasma tube center is 

= ^4- (37) 

Since Eq. (36) is difficult to use, let us construct its acceptable approximation. For this 
purpose, consider the behavior of the Lorentz factor 70 at low and high currents Jz- If the 
current Jz is low, then we can neglect the second term in Eq. (33) and obtain 

7o = lext, Jz 0. (38) 

As it must be, the screening of the external longitudinal electric field disappears at low 
currents. If, alternatively, the current Jz is high, then we can neglect the first term in 
Eq. (33): 

7o = lext(-r) , Jz > Jo- (39) 



■Jo 

Here, we introduced the screening current 

Jo = -lexU (40) 

a 

in exceeding which the screening of the longitudinal electric field becomes significant. If 
^ext ~ 10^ and ^0 ~ 0.1, then Jq ~ 10^. To construct the approximation, we will proceed 
from Eq. (39). Its shortcoming lies in the existence of a divergence when — > 0. To satisfy 
requirement (38), let us shift the argument in Eq. (39) by the screening current Jq: 

/ J\ 

7o^«7e.*(^l + 5^j . (41) 

It is easy to see that this approximation satisfies asymptotics (38) and (39). The accuracy 
of the approximation can be estimated by calculating (70 — 7o^)/7o, where 70 is the exact 
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solution (36). Obviously, this quantity depends only on the ratio Jz/ Jq- The dependence is 
plotted in the figure. We see that the approximation is rather accurate and the maximum 
deviation of 7q^ from 70 reaches 3.3% at J^/Jq ~ 1.62. This allows the approximation 

to be used instead of the exact solution 70 in the subsequent calculations. 

Using Eq. (37), we can derive the dependence of the longitudinal electric field on the 
longitudinal current J^: 

^(0) _ ^ext I ^ ^ 



J \ -4/3 



This dependence shows that the screening is switched on at currents Jz ^ Jo- Now, using 
Eq. (31), we obtain 

where we introduced the source 

equal to the effective source of particles when there is no screening of the external longitudinal 
electric field. Let ^ 10~^ and ^ext ~ 10^. The particle acceleration time to the stationary 
Lorentz factor is then Tgt = jext/Ew ~ 10^^, which gives an estimate of Qo ~ 10~^^ at ^0 ~ 
0.1. 



5. PARTICLE PRODUCTION BEFORE SCREENING 

Let us find the longitudinal distribution of the electron and positron densities and the 
corresponding total number of produced particles. For this purpose, we will take into account 
dependence (42) and rewrite Eq. (30) as 

dt dz 1 + Jz/Jo' ^ ' 

We will not write Eq. (29), which is an obvious corollary of the charge conservation law. 
The total source of particles, i.e., the total number of electrons and positrons produced per 
unit time in a plasma tube segment with a unit length, appears on the right-hand side of 
Eq. (44). At a low longitudinal current Jj, Eq. (44) transforms to 

dJz dNp 
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i.e., the electron-positron pair production rate is the same as that in an unscreened external 
electromagnetic field. If, alternatively, the current exceeds Jo, then 



(45) 



dt dz 

In this case, the pair production rate reaches some saturation value and the dependence on 
the total longitudinal current disappears. Hence, at a high longitudinal current Jz ^ Jo, 
when the longitudinal electric field screening processes come to the fore, the total electron- 
positron pair production rate coincides with that in the absence of screening but with a fixed 
longitudinal current J^ = Jq. Thus, here we observe the saturation of the particle source 
and the critical current Jq may be called the saturation one. 

By differentiating Eq. (44) with respect to time by taking into account the charge con- 
servation law (29), we can exclude the linear charge density from consideration and consider 
the equation for the longitudinal current J^: 

n 1 2Qo dJz 

" ^ (1 + jjj.y dt • 

For this equation, we should consider the Cauchy problem with the initial conditions 



J 



2lt=0 



= J 



init 



dt 



init 



Qj\jinit 



t=0 



1 + Jf'^/Jo dz 



where 



N: 



init 



•p "'Plt=0" 

The initial hnear charge density A^^"** and the initial current J^"** are some arbitrary func- 
tions of the longitudinal coordinate z. 

Let initially J*™* ^ Jq. As long as Jz < Jo, the plasma multiplication will be described 
by a linear hyperbolic equation: 

dJ. 



□ll Jz = -2(5o- 



dt ' 



Jz < ^0- 



Its solution is 



Jz{z,t) = 



z+t 



I loiQoVt' - {z - z'y) 



dJ., 



dt 



{z')-QoJz\t=o i^')) dz' 

=0 / 



+ 



z+t 

I 



z-t 



Z') 



= h {Qo^/t^ -{z- z'f) J.|,=o (^') dz' 



(46) 



+ Jz\t=0 {z-t)+ Jz\t^Q {z + t) 
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where loiz) and Ii{z) are, respectively, the zeroth- and first-order Bessel functions of an 
imaginary argument. Using solution (46) and the charge conservation law (29), we obtain 
the linear charge density 

t 

N,{z, t) = N^;^\z) - j ^^^^ dt'. (47) 



Hence we can derive the expressions for the linear electron and positron densities 

7 ± TV 

= (48) 
For the case under consideration, we should take 

jinit jyinit _j_ jyinit jyinit jyinit jyinit 

where N^^^ is defined by the formula 

Nr'-^^0{h-\z\). (49) 

Here, 9{x) is the theta function and li is the distance from the primary electron-positron pair 
produced by a photon from the external cosmic gamma-ray background at which secondary 
electron-positron pairs are produced [1]. 

It makes sense to consider solution (46) at currents that arc low compared to the 
screening current Jq. To estimate the time until which this can be done, note that, being 
symmetric relative to the initial point z — 0, solution (46) at long times has a maximum at 
this point. The required time that we will denote by tgcr can then be determined from the 
equation 

Jz{0,tscr) = Jo- (50) 

We will assume that i » ti Zi, where ti is the time corresponding to the distance h. The 
equality 

UO,t) = 3goe'^°*(/o(W + /i(Qot)), t > n, (51) 

then follows from Eq. (46). We will also assume that Qot ^ 1. Under the condition ti ~ Tgt, 
it follows from Eq. (31) that t ^ Ti^q- Since ^ 1, this condition does not restrict us in 
any way if f ^ ti. The asymptotics of the Bessel functions of an imaginary argument [32] 




z — >■ oo. 
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then allows Eq. (51) to be rewritten as 

Now, Eq. (50) can be brought to the equivalent form 



Cscr — ^scr ~l~ 2 ^^^scri (^^) 

where 

Acer = In — ^— • 
6Q0 

At Jo ~ 10^ and Qq ~ 10~^^, we have A^cr ~ 50. At A^^r iS> 1, Eq. (52) has an approximate 
solution, 

Cscr ~ ^scr ~^ '^^^ ^scr ■ 

For the source Qq and the linear longitudinal screening current density Jo mentioned above, 
we have ^scr ~ 52, which corresponds to Qotgcr ^ 1- The sought-for screening time is then 



tscr — „^ , (53) 



with tscr ~ 10^'' - 10^^ The ratio 



scr 



Tst 

is approximately equal to tscr/ Tst ~ 5 > 1. It remains to determine the total number A^. 
of electron-positron pairs produced before screening. Integrating the equation 

written for the total number of electron-positron pairs, 

oo 

N^{t)^ J N±iz,t)dz, (54) 

— oo 

and choosing A^|,"** = 3 as the initial condition (see Eqs. (49) and (54)), we obtain 

jyscr _ 2g2Qotscr 

Finally, using Eqs. (52) and (53), we can derive the final expression 
At Jo ~ 10^ Qo ~ 10-l^ and ~ 0.1, we obtain A^f ~ lO^^. 
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6. PARTICLE PRODUCTION DURING SCREENING 



After the time Tgcr, the total source of electron-positron pairs becomes constant due to 
the saturation effect (see Eq. (45)). However, we should take into account the fact that 
this source is constant only on distance scales that are small compared to the scales of the 
change in external magnetospheric electromagnetic field coincident with the neutron star 
radius. In the preceding section, this did not require a special allowance, because tscr ^ Rs- 
In contrast, the plasma generation processes will subsequently develop on scales of the order 
of Rs — the scales that the plasma tube length will eventually reach. Let us introduce the 
time 



after which (from the appearance of a primary electron-positron pair) the plasma generation 
will occur under saturation. It can be assumed that to ~ tscr- Let us designate the total 
source of electron-positron pairs as 



We will assume it to be some arbitrary function of the longitudinal coordinate z. This can 
be done, because Qsiz) depends only on the external parameters of the magnetospheric elec- 
tromagnetic field but not on the electron and positron densities in the plasma tube. Recall 
that the production point of a primary electron-positron pair fixes some magnetic field line. 
The produced secondary particles will move almost along the field line under consideration 
and their positions can be specified by some longitudinal coordinate z. Fixing the field line 
for known external electromagnetic fields completely determines the function Qs{z). Let 
us also introduce the longitudinal coordinate Zq corresponding to the production point of a 
primary electron-positron pair. Without loss of generality, we will set Zq = 0. To calculate 
the linear electron and positron densities, it is more convenient to directly use Eq. (28) and 
to write 



In this case, we consider to to be the initial time to be able to write the initial conditions 
derived from the calculations of the preceding section as the conditions at t — 0. Direct 



to = n + 1. 



sen 



Qs{z) = QqJq. 




(56) 
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integration of Eqs. (56) gives the linear densities of positrons 



z z 

9{to + t+)J + 9{to + t-)J 



+ e{to-t^)J - J 



and electrons 



N_{z,t)^N_\^^,{U) + 



(to+t-)/2 -t- 
(to+t-)/2 z 

QMdy 

t- -t- 



t+ (to+t+)/2 



t+ t+_ 



QMdy 



{to+t+)/2 



and the longitudinal current 



Jzi 


^,^) = ^+Uo(-^-)- 






{to+t+)/2 


t+ 


+ 




-t+)j + 9{to 




-t- 


(to+t+)/2 



9{to + 



-(to+t-)/2 
-(to+t-)/2 t+ 



(57) 



Qii{y)dy. 



We will not need the expression for the linear charge density Np. The expressions for and 
derived from Eqs. (46)-(48) after the substitution t — tgcr should be taken as -A^+|t=o (-^) 
and Ar_|^^o (z). 

The total number of produced electron-positron pairs is equal to half the value obtained 
by integrating Eq. (57) over the longitudinal coordinate z: 



where 



jysat ^ 



to+t —to to+t 

tj + j{h + y) + j{h-y) 



QMdy 



(58) 



—to—t —to—t to 

is the total number of electron-positron pairs produced after the onset of longitudinal electric 
field screening. In this formula, the number of pairs produced before the onset of screening 



19 



acts the number of electron-positron pairs at the initial time, i.e., -/Vs|t=o ~ ^^^^ (^®® 
Eq. (55)). It remains to use Eqs. (40) and (43) and to write the function 



where jext and Tgt are considered as functions of the longitudinal coordinate z. Recalling 
the dependence of the particle acceleration time on the stationary Lorentz factor and the 
external longitudinal electric field -Bl^*, we ultimately find 



Tpext 

Qe(^) = (59) 



At * - 10-^ we have Qj:{z) - 10" 



7. CALCULATION FOR A DIPOLE MAGNETIC FIELD 

The subsequent calculation of the specific number of electron-positron pairs requires 
knowing the dependence of the external longitudinal electric field on the parameter z 
along a magnetic field line. This dependence is completely determined by the production 
point of a primary electron-positron pair in a neutron star magnetosphere. However, the 
following considerations can be used for qualitative estimates. Let the primary electron- 
positron pair be produced at a distance r from the stellar center. The total length of the 
field line fixed by the pair production point will then also be of the order of r. In this case, 
the length of its part located above the neutron star surface is equal in order of magnitude 
to r — Rs- For simplicity, the produced secondary particles can be assumed to move along 
the radial coordinate. We see from the expressions for Deutsch's electromagnetic field [28] 
that the electric field changes along the radial coordinate as a power law. For the subsequent 
estimates, we will assume that 

E^^\r) ~ (60) 

where m is the magnitude of the neutron star magnetic moment, kn = 1/Rl is the wave 
number corresponding to the angular frequency Q of neutron star rotation, and Rl = c/Q 
is the light cylinder radius written in dimensionless units. 

Let us first calculate the total number of produced electron-positron pairs for the general 
model power dependence 

Qj:{z) = Ao{z + Zo)-''. (61) 
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We used the fact that under the assumptions described above, r = z + zq, where zq is 
the distance from the neutron star center at which the primary electron-positron pair is 
produced from which the plasma tube is subsequently formed. We will also assume that 
u > 2. The substitution of dependence (61) into Eq. (58) and subsequent integration then 
yield the expression 



jysat ^ 



{u-l){u-2) 



(;^o-io-i)~"+' + (^o + *o + ^)""+' 

(62) 

■(^0-^0)-"+' -(^0 + ^0)-"+' 



We will use the fact that the switch-on time of longitudinal electric field screening (in 
dimensionless units) is small compared to the neutron star radius. Consider the case where 
the distance from the neutron star center at which the primary electron-positron pair is 
produced is large compared to the stellar radius. We can then write the inequalities 

If, in this case, we take 

t = Zq — to — Rs, 

then it follows from Eq. (62) that 

= ^ ^i?5"+'- 

^ {u-l){u-2) ^ 

The latter quantity is the number of electron-positron pairs produced in the regime of 

saturation by the instant the plasma tube reaches the neutron star surface. Let us now 
return to Eq. (59). If we take into account Eq. (60), then = kamRg /2a and u = i 
correspond to it. Then, 

^ 12 q; 

is the characteristic number of electron-positron pairs produced in the plasma tube after the 
onset of longitudinal electric field screening, so that the plasma is generated in the regime 
of saturation of the total source. It is easy to see that Njf^^ depends only on the dipole 
magnetic moment of the neutron star and on the angular frequency of its rotation. If the 
dipole magnetic moment m of the neutron star is expressed in terms of its radius Rs and 
the magnetic field strength Bq at its magnetic pole, then 
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For Rs ~ 2.6 x 10^^ (10 km in dimensional units), Rl — ^O^Rs, and Bo ~ 0.01, we have 
7V|,"* ~ 10^^. If, alternatively, the magnetic field strength at the magnetic equator of the 
neutron star is meant by Bq, then the right-hand side of Eq. (63) should formally be doubled. 

Thus, we can write the total number of electron-positron pairs produced in the plasma 
tube from the instant of primary pair production by a photon from the external cosmic 
gamma-ray background to the instant the plasma tube reaches the neutron star surface, 

where Nj^^ and N§^^ are defined by Eqs. (55) and (63), respectively. Given the inequality 

we can ultimately write 

It remains to verify that the effective source Q^^^ can be used over the entire plasma 
generation time. This requires that the condition A^^ ^ 1, where r = Tst is the delay time, 
be met. This inequality holds if the particle Lorentz factor exceeds 10^ [1]. The longitudinal 
current can be estimated from the formula 

jvrsat 

J. < 2-^. (64) 

In order of magnitude, ~ 10^^ — 10^^, which corresponds to Jz/Jo ~ 10^ — 10^. Conse- 
quently, the longitudinal electric field decreases by three orders of magnitude compared 
to the initial external electric field, which corresponds to ~ 10~^ for E^^^ ~ 10~^. It 
follows from Eq. (41) that the Lorentz factor decreases by only one order of magnitude 
compared to the initial one, 7q ~ Q.l'fext, and is 7q ~ 10^ for ^ye^t ~ 10^. 
Note the interesting relation 

iVv 1 



ugjVs 16' 

where Vs = 47ri?|/3 is the volume of the neutron star and Ugj is the Goldreich- Julian density 
near its surface. It shows that the initial filling of the magnetosphere is very efficient. Only 
10—100 absorbed photons from the external cosmic gamma-ray background are enough to fill 
the inner neutron star magnetosphere. The lightnings forming in this case contain a number 
of particles comparable to that in the inner magnetospheric regions completely filled with 
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plasma. If, alternatively, we consider the entire magnetosphere up to the light cylinder Rl, 
then the number of particles in it is equal in order of magnitude to 3nGjVsln{RL/ Rs). 
At Rl/Rs ~ 10^ — 10^, about 300 — 500 absorbed photons are needed to fill the entire 
magnetosphere. 



8. THE PLASMA TUBE RADIUS AND MULTIPLICITY 

Let us estimate the multiplicity 

Tl 

Am = , (65) 

where n is the number density of electron-positron pairs in the plasma tube and ncj is the 
Goldreich- Julian density equal in dimensionless units to 

\Q-B\ 

riGj = 



27ra 

For our estimates, we can neglect the dependence on the angle between the vectors fl and 
B present in this formula and write 

''^'^2^R-, ^''^ 
We know typical values of the longitudinal current J^. To estimate the number density n, 
the radius of the forming tube should be known. 

It is easy to estimate the radius of the plasma tube if the characteristic velocity of its 
boundary is known. This velocity is determined by two factors. The first factor is the 
electric drift of particles at the tube boundary in crossed electric and magnetic fields (see 
Eq. (1)). In order of magnitude, the electric drift velocity can reach E/B ~ 10~^ (in units 
of the speed of light). The second factor is the motion of the tube boundary due to the finite 
photon mean free path If. The point is that the photon emitted by a particle almost along 
the tangent to the magnetic field line produces an electron-positron pair on a different field 
line shifted relative to the initial field fine by the distance 

due to the curvature of the field line itself. This distance corresponds to the orthogonal 
velocity of the boundary 

'/ 
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For our estimates, as Z/ we can take l^ax — PXmax, where Xmax is related to the minimum 
longitudinal Lorentz factor of the produced particles as 7inin = 1/Xmax [17]. We will then 
obtain the following qualitative estimate: 

V, = (67) 

/mm 

For 7inin ~ 10^, we obtain v± ~ 10~^ from this estimate. We see that the finite mean free 
path can increase considerably the velocity of the boundary compared to what the ordinary 
electric drift velocity can give. The time dependence of the tube radius will be given by the 
formula 

R = v±t, 

where, for simplicity, the velocity of the boundary is assumed to be constant. For this 
estimate, the tube radius on time scales t ~ Rs will reach 

R < 100 m. (68) 

To estimate the multiplicity, wc will now take the electric current estimate from Eq. (64), 
in which a tube length equal to Rs was taken. The corresponding time is t = Rsf^. The 
number density of electron-positron pairs is then 

2 J. 



n 



vl ttRI' 



Using Eqs. (63)-(66), we ultimately obtain 

A. 



3vl 



Since we took Bq as the magnetic field on the magnetic equator of the neutron star, Eq. (63) 
was modified in accordance with the remark made after it. We see that to determine typical 
values of the multiplicity, it will suffice to know only the expansion velocity of the forming 
plasma tube. The lower multiplicity limit can be derived using Eq. (67) after the substitution 
7min ~ 100. The upper limit can be obtained if a typical drift velocity E/B ~ 10^"^ is taken 
as the velocity of the boundary. Thus, we have 

10^ < A™ < 10^ 



Obviously, the upper limit is grossly overestimated, because the finiteness of the photon 
mean free path is disregarded, and the actual multiplicity is closer to the lower limit. 
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9. CONCLUSIONS 



The absorption of a high-energy photon from the external cosmic gamma-ray background 
in the inner neutron star magnetosphere triggers the generation of a secondary electron- 
positron plasma and gives rise to a lightning — a lengthening and simultaneously expanding 
plasma tube. It propagates along magnetic field lines with a velocity close to the speed of 
light. The orthogonal velocity of the tube side boundary can reach 0.01 of the speed of light. 
Such a velocity is provided by the pair production on magnetic field lines shifted relative to 
those along which the photon-emitting particles move. 

The shift is determined by the finite photon mean free path and the curvature of the 
magnetic field lines. The existence of a strong longitudinal electric field provides the accel- 
eration of particles and the emission of curvature and synchrotron photons by them, which, 
in turn, also produce electron-positron pairs. Initially, the plasma is produced in the ex- 
ternal vacuum electromagnetic field, while there is no screening of this field. After a time 
tscr ~ 10~^ — 10~^ s, the electric current in the lightning becomes equal to the screening cur- 
rent Jo ~ 10^^ A. Note that the screening current is comparable to the characteristic electric 
current fiowing through the polar region of the magnetosphere of an ordinary radio pulsar. 
About 10^^ electron-positron pairs are produced by the instant the screening current in the 
lightning is reached. Subsequently, the dynamical screening of the longitudinal electric field 
through the growing electron-positron generation rate begins. It is provided not by static 
charge separation but by electric current growth in the lightning. We emphasize that the 
production of an electron-positron plasma in a lightning is essentially nonstationary, while 
the dynamical electric field screening in the stationary case is impossible. 

The plasma is mainly produced in the regime of screening, when the longitudinal electric 
field in the lightning and the source of electron-positron pairs are determined self-consistently 
from the following condition: the longitudinal electric field depends on the electron-positron 
pair generation rate, which itself depends on the electric field in the lighting. The time it 
takes for the lightning to reach the neutron star surface may be considered as its lifetime. 
Subsequently, the screening mechanism can change due to the infiuence of the neutron star 
surface. The plasma production can even cease, for example, due to the onset of static 
electric field screening. The electric current fiowing in the lightning when it reaches the 
neutron star surface is higher than the screening current by two or three orders of magnitude. 
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Although the longitudinal electric field by the end of lightning formation decreases by three 
orders of magnitude compared to the initial vacuum electric field, the particle Lorentz factor 
decreases by only an order of magnitude and reaches 10^—10^. This value is much larger than 
the characteristic Lorentz factor of the secondary-plasma particles in the magnetosphere of 
an ordinary pulsar. 

If the production of an electron-positron plasma is assumed to be efficient in a magnetic 
field higher than 10^ G, then the distance of (10 — 100)Rs in exceeding which the plasma 
generation begins to be suppressed corresponds to a neutron star surface magnetic field 
of 10^^ — lO-*^^ G. Consequently, the lightning length can reach 1000 km. The lightning 
radius is 100 m. Interestingly, the lightning radius is comparable to the radius of the polar 
regions in a neutron star magnetosphere formed by open field lines. The number of electron- 
positron pairs produced in the lightning in its lifetime reaches 10^^. This corresponds to the 
constraints on the multiplicity > 10^, which is comparable to or even larger than that 
in the polar regions of ordinary pulsars. 

We see that the parameters of the generated plasma are comparable to those of the 
electron-positron plasma in the polar regions of ordinary radio pulsars. Indeed, the multi- 
plicity is comparable to that in the magnetosphere of a radio pulsar. The particle Lorentz 
factor exceeds considerably that of the secondary plasma in pulsars and magnetars. The 
curvature of the magnetic field lines is comparable to that in the polar cap of a radio pulsar 
if the lightning is formed in the region of open field lines or exceeds it considerably if the 
lightning is formed in the region of the magnetic equator. In addition, the lightning width is 
comparable to the polar cap width of a radio pulsar. In view of the existing analogy, these 
factors suggest that the radio emission can be generated by individual lightnings. 

If we consider a neutron star located at a distance of the order of the distance to a 
typical pulsar, then we can assume the possibility of observing radio bursts from lightnings at 
existing radio astronomy observatories. Moreover, one may expect the radio fiux density in a 
lightning to exceed the flux density from ordinary radio pulsars due to larger particle Lorentz 
factor, multiplicity, and plasma density. One of the factors determining such parameters is 
the following difference between the plasma generation mechanisms in a lightning and in 
the polar cap of a radio pulsar. If we consider an ordinary radio pulsar with a relatively 
weak (compared to the critical one) magnetic field, then, apart from the primary particles 
accelerated in the polar gap, there are at least two more generations of secondary particles: 
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the first generation of secondary particles produced by the curvature photons from primary 
particles and the second generation of particles produced by the synchrotron photons emitted 
through the transition of first-generation particles to the zeroth Landau level. However, a 
nearly vacuum electric field is available only in the polar gap up to distances of about 100 m 
above the neutron star surface. Above the gap, the electric field is screened. In contrast, 
in the case of a lightning, even despite the dynamical screening, we have a fairly significant 
longitudinal electric field that continuously accelerates the newly produced particles. These 
particles, whatever generation they formally belong to, immediately begin to play the role of 
primary particles after their acceleration. Obviously, the electric field in a lightning cannot 
strictly become zero. If this occurred, then the electric field screening would immediately 
disappear and the field itself would become close to Deutsch's vacuum field. Consequently, 
the electric field is always nonzero and the plasma is generated continuously as long as the 
lightning exists. 

Stationary electron-positron plasma outfiows from the magnetosphere are responsible for 
the radio emission from pulsars, but the plasma generation in a lightning is nonstationary. 
Consequently, the radio emission from lightnings, if it can be observed, must be connected 
with the observational manifestations of nonstationary radio sources associated with neutron 
stars. The formation time of radio emission is limited by the lightning lifetime. If the width 
of a radio burst is assumed to be comparable in order of magnitude to the ratio of the 
lightning length to the speed of light, then we have a burst width of 3.3 ms for a lightning 
length of 1000 km. Radio bursts of comparable width are observed from rotating radio 
transients. For example, the duration of a single burst from RRAT J1819-1458 is 3 ms in 
observations at 1.4 GHz [3]. Here, we disregard the pulse broadening due to scattering, 
which is important in RRAT observations at lower radio frequencies [33]. As regards the 
fiux density, it reaches 10 Jy at 1.4 GHz [34], exceeding considerably that for typical radio 
pulsars. 

Thus, we can assume that the radio bursts from lightnings can be observed as those from 
RRATs. Investigation of the observational manifestations of lightnings in a neutron star 
magnetosphere and their possible connection with radio bursts from RRATs is beyond the 
scope of this paper and requires a separate consideration. 
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FIG. 1. Accuracy of the approximation of the Lorcntz factor 70 versus ratio of the electric 
current flowing in the lightning to the screening current Jq. 
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